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SECTION  I 


INTRODUCTION 

Reentry  vehicles  may,  depending  on  weather  conditions,  encounter  water  droplets 
and/or  ice  crystals  at  altitudes  below  50,000  ft.  The  concentration  of  condensed 
phase  water  (hydrometeors)  is  typically  on  the  order  of  hundredths  to  tenths  of  a 
gram  per  cubic  meter  for  ice  crystals  in  cloud  formations  and  can  be  on  the  order  of 
several  grams  per  cubic  meter  for  water  droplets  in  rain  fields.  Due  to  high  velo- 
cities a reentry  vehicle  can  experience  severe  erosion  damage  arising  from  encounters 
with  hydrometeors.  In  order  to  accurately  predict  erosion  damage  experienced  by 
specific  materials  it  is  necessary  to  know  hydrometeor  impact  mass  and  velocity.  The 
values  of  the  aforementioned  quantities  will  be  affected  by  phenomenon  associated 
with  interaction  between  the  shocked  gas  layer  surrounding  the  vehicle  and  the  hydrom- 
eteors traversing  this  layer;  these  effects  decrease  with  increasing  hydrometeor  size 
and  decreasing  shock  strength.  Consequently,  in  order  to  accurately  assess  erosion 
damage  caused  by  variable  size  hydrometeors  and  variable  vehicle  geometries  and  flight 
conditions,  it  Is  necessary  to  understand  and  be  able  to  quantify  hydrometeor/shock 
interaction  phenomenon.  Thus  an  experimental  program,  utilizing  a shock  tube  facility 
was  Initiated,  directed  at  determining  the  effects  of  shocked  gases  on  hydrometeors. 

To  date  this  program  has  been  successful  in  obtaining  data  relative  to  water  droplet/ 
shock  interactions  and  there  is  an  on-going  effort  directed  at  establishing  a data 
base  relative  to  Ice  crystal/shock  Interactions.  The  present  report  documents  the 
results  and  data  correlations  obtained  from  the  water  droplet  tests. 

There  have  been  a number  of  previous  investigations  analytical  and  experimen- 
tal In  nature  on  particle  gas  Interactions  that,  because  of  their  relevance,  should 
be  cited.  Early  studies  have  been  concerned  with  drag,  heat  transfer,  mass  transfer, 
and  break-up  of  particles  in  moderate  to  high  energy  gaseous  flow  fields.  Analyses 
of  particle  heating  and  drag  for  rocket  nozzle  environments  are  presented  In  Refer- 
ences 1 through  3.  References  4 and  5 contain  treatments  of  mass  transfer  effects 
arising  from  chemical  reactions  for  combustion  chamber  as  well  as  nozzle  environ- 
ments. In  the  aforementioned  references  drag  effects  are  modelled  through  the  appli- 
cation of  drag  coefflcient-Reynolds  number  correlations  and  heat  transfer  is  modelled 
through  the  application  of  Nusselt  No.  correlations  with  Reynolds  No.  and  Prandtl  No. 
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More  recently  a number  of  analytic  and  experimental  Investigations  have  been 
conducted  on  droplet  and  Ice  crystal  phenomenon  behind  shock  waves  (see  References 
6 through  17).  References  6 through  8 treat  weak  shock  wave  effects;  references  9 
through  17  are  specifically  concerned  with  droplet  and  ice  crystal  shock  interactions 
for  reentry  flight  conditions.  The  earlier  references  (9  through  12)  present  experi- 
mental results,  data  correlations,  and  a limited  amount  of  analysis.  While  the  work 
described  in  References  13  through  17  is  purely  of  an  analytic  nature. 

The  results  documented  in  Reference  9 were  for  shock  Mach  numbers  of  3,  6,  9, 
and  11  at  pressures  of  1 atm,  0.46  atm,  and  0.184  atm,  the  lower  pressures  corre- 
sponding to  higher  Mach  numbers.  Water  drop  diameters  were  varied  from  500p  to  2500w, 
a range  that  is  representative  of  naturally  occurring  rain  drops.  Droplet  break  up 
times  and  droplet  displacement  histories  were  obtained  by  piecing  together  informa- 
tion from  single  frame  photographic  data.  In  addition  X-ray  photography  was  employed 
to  obtain  mass  histories.  It  was  found  that  the  displacement  data  could  be  corre- 
lated in  terms  of  the  parameters  X and  T where 

x/do 

(t  U2/d0)^7p[ 

droplet  displacement  from  unshocked  location 
gas  velocity  behind  shock 


X 

T 

x 

U, 


dQ  ■ unshocked  droplet  diameter 
p2  * gas  density  behind  shock 
* droplet  density 

and  the  correlating  function  is  given  by 
X * 0.8T2 

In  addition  It  was  found  that  at  a sufficient  time  Interval  after  passage  of  the 
shock  the  droplets  broke  up  catastrophically,  and  a dimensionless  break  up  time,  Tb, 
could  be  obtained  from  the  following  Weber  number  correlation 


Tfe  « 45/We1/" 


where 


We 


°2U*dc 


surface  tension  evaluated  at  conditions  upstream  of  shock 


Prior  to  times  corresponding  to  Tb>  mass  loss  was  observed  due  to  stripping.  Appli- 
cation of  the  results  of  Reference  10  Is  restricted  due  to  the  limited  range  of  test 
conditions.  In  that  study  only  a single  sized  droplet  was  tested  (approximately  500p), 
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for  the  Mach  No.  range  3-4.3;  upstream  pressures  ranged  from  1 to  2 atm.  At  a Mach 
No.  of  9 the  free  stream  pressure  was  0.014  atm;  the  altitude  corresponding  to  this 
pressure  is  well  above  that  for  which  atmospheric  water  droplets  will  exist.  The 
results  for  drop  displacement.  Reference  10  with  the  exception  of  the  high  Mach  No. 
low  pressure  case,  are  in  reasonable  agreement  with  those  in  Reference  9.  For  the 
high  Mach  No.  case  the  droplet  acceleration  rates  were  higher;  however,  as  was  men- 
tioned the  test  conditions  for  this  case  are  not  representative  of  reentry  flight 
conditions. 

A method  for  predicting  water  droplet  trajectories  through  reentry  vehicle 
shock  layers  including  mass  loss  effects  due  to  evaporation  is  presented  in  Reference 
11;  particle  breakup  phenomena  is  also  discussed  in  this  reference.  The  analysis 
utilizing  numerous  simplifying  assumptions,  starting  with  the  equations  of  motion, 
develops  closed  form  solutions  for  droplet  coordinates  (x,y)  as  a function  of  time 
as  the  droplet  traverses  the  shocked  gas  layer.  As  a result  of  the  simplifying  assump 
tion  the  equations  of  motion  are  uncoupled  from  the  mass  loss  equation  which  is 
solved  numerically.  In  References  13  and  14  an  improved  model  for  heat  transfer  is 
used  to  solve  for  the  mass  (size)  variation  with  time  and  the  equations  of  motion  are 
coupled  to  the  mass  loss  rate  equation.  It  is  pointed  out  that  the  solutions  in 
Reference  12  are  more  general  than  those  in  References  13  and  14  from  the  standpoint 
that  the  former  treats  regions  of  stagnation  and  conical  flow  whereas,  strictly  speak- 
ing, the  results  in  References  13  and  14  are  valid  only  along  the  stagnation  stream- 
line. 

An  extension  of  the  data  base  and  analysis  reported  in  Reference  9 is  given 
in  Reference  11;  essentially  Reference  11  presents  results  and  data  correlations 
covering  smaller  size  droplets  than  the  lower  limit  for  tests  reported  in  Reference 
9.  The  principle  correlating  parameter  in  the  aforementioned  References  is  Weber  No. 
and  Reference  11  contains  correlations,  including  results  of  other  investigations,  of 
droplet  displacement  and  break  up  data  over  a five  order  of  magnitude  range  of  Weber 
No.  The  large  scatter  about  the  suggested  correlations  indicates  that  possibly  other 
parameters,  in  addition  to  Weber  No.,  are  of  importance. 

Reference  15  represents  the  first  attempt  to  model  hydrometeor/shock  inter- 
action phenomena  by  solving  the  governing  flow  equations  on  a local  level  relative 
to  the  droplet.  A two  layer  boundary- layer  model  is  employed  to  describe  heat  and 
mass  transfer  events  at  the  surface.  The  results  indicate  that  heat  and  mass  trans- 
fer effects  are  less  than  those  predicted  previously  by  applying  global  Nusselt  No. 
correlations  (i.e.,  References  13  and  14).  It  Is  felt  that  the  results  obtained 
using  the  more  fundamental  approach  (I.e.,  Reference  15)  showing  smaller  shock  layer 
effects,  are  caused  at  least  In  part  by  neglect  of  the  blocking  effect  on  heat  and 
mass  transfer  In  Reference  13  and  14. 
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A sophisticated  procedure  for  predicting  particle  trajectories  including  effects 
of  heat  transfer,  mass  transfer,  and  particle  breakup  is  presented  in  Reference  16. 

The  results  in  this  reference  are  based  on  a programmed  method  of  numerically  solving 
the  governing  equations  without  relying  on  simplifying  assumptions.  Drag  is  treated 
by  employing  Reynolds  No.  correlations  and  heat  transfer  is  treated  in  a global  man- 
ner by  applying  correlations  for  Nusselt  No.  which  include  mass  transfer  effects.  Sub- 
sequently this  method  has  been  applied  to  treating  the  effects  of  shock  layers  on 
hypervelocity  impact  induced  erosion  (Reference  17),  coupled  with  thermal  ablation. 

The  present  investigation  has  been  experimental  and  represents  an  extension 
of  the  work  described  in  Reference  12.  The  time  duration  from  droplet/shock  encounter 
and  test  conditions  for  the  present  study  were  chosen  to  give  results  pertinent  to 
the  stagnation  region  of  reentry  vehicles.  The  bulk  of  the  data  was  obtained  utiliz- 
ing an  extremely  high  speed  multi-framing  camera  in  order  to  obtain  a high  resolu- 
tion of  events  in  the  immediate  region  behind  a shock  front.  The  overall  program 
objectives  are  to:  (1)  determine  reentry  vehicle  shock  layer  effects  on  hydrometeor 

breakup  and  (2)  determine  reentry  vehicle  shock  layer  effects  on  hydrometeor  impact 
velocity.  The  approach  adopted  to  accomplish  the  aforementioned  objectives  consists 
of  (1)  utilization  of  a high  energy  (NASA  AMES  EAST  Facility)  shock  tube  to  generate 
test  conditions,  (2)  coupling  droplet  generator  to  shock  tube  test  section,  (3)  uti- 
lization of  high  speed  multiframing  (Hadland  Imacon)  camera  to  obtain  data  - 8 to  16 
total  frames  at  framing  rates  up  to  20  x 10e  per  second. 

It  is  pointed  out  that  there  are  two  advantages  of  a multiple  framing  capa- 
bility relative  to  that  of  a single  framing  capability,  namely:  (1)  greater  effi- 
ciency in  data  acquisition,  and  (2)  improved  accuracy.  The  first  advantage  stems 
from  the  fact  that  multiple  data  points,  corresponding  to  different  distances  between 
shock  front  and  hydrometeor  can  be  obtained  for  each  test  condition,  with  a single 
frame  camera  only  at  a rate  of  one  point  per  shot.  The  second  advantage  stems  from 
the  fact  that  error  sources  exist  in  a data  set  obtained  with  a single  frame  camera 
that  are  automatically  eliminated  with  a multi-frame  camera;  examples  are  (1)  uncer- 
tainty In  Initial  droplet  position  - resulting  in  errors  of  measured  droplet  dis- 
placement, (2)  uncertainty  In  Initial  droplet  size  - resulting  in  erroneous  scaling 
utilized  In  dimensionless  correlations,  (3)  uncertainty  in  shock  location.  That  is. 

In  order  to  deduce  a sequence  of  events  experienced  by  a single  droplet  from  single 
frame  data  It  Is  assumed  that  variation  in  test  conditions,  droplet  size,  and  drop- 
let location  are  negligible  from  test  to  test.  Furthermore,  it  implies  an  extremely 
high  degree  of  accuracy  of  Instruments  employed  to  locate  shock. 

The  following  sections  describe  the  test  equipment,  test  results,  data  cor- 
relations, and  the  application  of  the  correlations  to  the  prediction  of  shock  layer 
effects  on  erosion  attenuation. 
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SECTION  2 


EXPERIMENTAL  APPARATUS 

The  essential  components  of  the  experimental  apparatus  are  the  shock  tube  and 
driver,  the  water  droplet  generator,  high  speed  camera,  light  source  and  accompanying 
optics.  Mounted  within  the  shock  tube  is  a test  section  with  windows  and  fixtures 
needed  to  accommodate  the  light  source,  optics,  and  droplet  generator.  The  following 
sections  describe  the  aforementioned  components. 

2.1  SHOCK  TUBE,  DRIVER,  AND  TEST  SECTION 

The  NASA  AMES  Electric  Arc  Shock  Tube  (EAST)  Facility  was  employed  in  this 
study  to  generate  test  conditions.  A drawing  of  the  facility  and  its  peripheral  ele- 
ments is  shown  in  Figure  2-1.  The  shock  tube  has  a 4 inch  internal  diameter  and  an 
8 Inch  outside  diameter.  The  outer  cross  section  of  the  test  section  is  a square 
having  side  length  seven  and  one-half  inches.  The  inner  cross  section  is  circular 
and  matches  that  of  the  shock  tube.  The  overall  length  of  the  tube  is  forty  feet. 

The  downstream  end  discharges  into  a 296  ft3  dump  tank.  The  tube  contains  ports  for 
ion  probes,  pressure  gauges  and  photomultipliers.  The  positions  of  these  ports  are 
given  in  Figure  2-1.  The  dimension  given  for  each  station  is  the  distance  from  that 
station  to  the  driver  exit  plane.  A membrane  which  ruptures  during  discharge  is 
located  at  the  driver  exit  plane. 

The  ion  probes,  designed  and  made  by  Ames  Research  Center,  measure  shock  posi- 
tion as  a function  of  time.  The  sensor  is  a twin  electrode  probe  mounted  in  the  wall 
of  the  tube  with  a small  dc  voltage  impressed  across  the  electrodes.  A power  supply, 
amplifier,  and  gate  pulse  generator  complete  the  instrument.  For  additional  details 
on  this  instrument  see  Reference  18.  The  pre-shot  shock  tube  pressure  is  obtained 
from  a Wallace  and  Tlernan  gauge  located  at  F-|.  A Kistler  gauge  for  obtaining  pres- 
sure during  a shot  Is  located  at  L-3. 

The  energy  producing  the  shock  is  produced  by  an  electric  discharge  In  the 
driver.  The  driver  chamber  Is  conical  with  the  high  voltage  electrode  located  in  the 
center  of  the  small  end  of  the  cone;  a ring  shaped  wall  flush  electrode  is  located  at 
the  larger,  downstream  end  of  the  driver  chamber.  Energy  for  the  driver  is  supplied 


2-1 


by  a 1250-p  Farad  capacitor  bank  rated  at  106  joules  when  charged  to  40  Kilovolts. 
Helium,  hydrogen  and  argon  driver  gases  have  been  employed;  in  the  present  tests 
hydrogen  was  used.  It  should  be  pointed  out  that  the  basic  shock  tube  independent 
variables  for  a given  test  are  capacitor  charging  voltage  and  shock  tube  pressure; 
the  former  quantity  determines  the  capacitor  storage  energy.  An  oscilloscope  trace 
of  the  current- voltage  history,  during  discharge  for  a typical  test  is  shown  in 
Figure  2-2.  Figure  2-3  shows  the  facility  operating  limits  during  the  present  pro- 
gram. The  curve  represents  an  upper  limit  on  attainable  velocity  for  a given  pres- 
sure; the  capacitor  voltage  drop  for  this  curve  is  nominally  15  - 17.5  KV.  Although 
the  capacitor  is  designed  to  higher  voltage  drops;  higher  voltages  do  not  produce 
higher  gas  velocities  due  to  increased  losses.  Figure  2-3  also  shows  the  facility 
operating  conditions  that  are  currently  attainable  using  a modified  driver.  Further 
details  on  the  NASA  Ames  EAST  Facility  are  given  in  References  19  and  20. 

2.2  DROPLET  GENERATOR 

Consistent  with  test  objectives  of  the  present  program  a system  was  developed 
capable  of  generating  water  droplets  ranging  in  sizes  down  to  50  pm  diameter  spheres. 
Accurate  positioning  as  well  as  size  control  was  a requirement.  Depth  of  field  and 
field  of  view  constraints  necessitated  that  the  droplet  trajectories  fall  within  a 
predetermined  location  not  greater  than  0.1  mn  wide  and  0.1  inn  deep. 

Two  droplet  generating  techniques  were  considered.  The  first  was  originally 
developed  by  Abbott  and  Cannon  (Reference  21)  at  the  National  Center  for  Atmospheric 
Research.  The  Abbott  and  Cannon  generator  utilizes  a wire  which  Is  repeatedly  In- 
serted into  and  withdrawn  from  a water  reservoir.  During  each  withdrawal  a cylin- 
drical filament  of  water  is  removed  from  the  surface.  When  the  fllamenter  reaches  a 
certain  length  it  breaks  off  and  contracts  into  a spherical  drop.  The  size  of  the 
droplet  is  a function  of  the  diameter  of  the  wire,  and  angle  between  the  wire  and 
surface  of  the  reservoir,  and  the  depth  of  penetration  of  the  wire  into  the  reservoir. 
The  droplet  once  formed  is  acted  on  by  gravity.  This  method  Is  extremely  versatile 
and  can  produce  droplets  ranging  In  diameter  from  several  hundred  microns  to  sub- 
microns. A disadvantage,  however.  Is  that  due  to  slow  settling  velocity  (0.2  ft/sec 
for  50  pm  diameter  droplets),  stray  air  currents  will  cause  large  amounts  of  drift; 
this  greatly  reduces  the  probability  of  obtaining  droplets  at  the  desired  location. 

For  this  reason  a pulsed  jet  type  of  generator  was  used  (see  References  22  and  23). 

In  these  generators  water  In  a pressured  reservoir,  with  a feed  line,  exits  a small 
diameter  nozzle  or  orifice.  Pressure  pulses  are  obtained  through  an  amplifier- 
oscillator  system  which  feeds  a signal  to  a piezoelectric  crystal  either  In  contact 
with  the  reservoir  surface  or  the  liquid  Itself.  The  pulses  Introduces  an  Instability 
In  the  free  jet  exiting  the  reservoir  which  results  In  breakup  and  subsequent  droplet 
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formation.  This  type  of  generator  can  generate  droplets  at  a greater  rate  (TO3  - lO” 
droplets/sec)  than  the  Abbott  generator  and  at  greater  velocities  (approximately  50 
ft/sec  for  50  umeter  droplets).  The  size  range  of  Interest  and  desired  accuracy  was 
achieved  with  this  type  of  generator. 

A schematic  of  the  droplet  generator  system  Is  shown  In  Figure  2-4.  Figure 
2-5  Is  a detailed  drawing  of  the  nozzle  assembly;  the  piezoelectric  crystal  Is  bonded 
to  the  upper  surface  of  the  nozzle  assembly.  A Heath  signal  generator  and  McIntosh 
amplifier  operating  off  of  line  voltage  activates  the  crystal.  The  Insulation  sec- 
tion and  components  attached  to  It  (see  Figure  2-4)  are  housed  In  a pressure  sealed 
container  which  Is  directly  connected  to  the  shock  tube.  Other  components  of  the 
system  are  a pressure  supply,  water  tank,  flow  control  valves,  rotameters,  and  fil- 
ters. Three  nozzles  having  diameters  of  0.0011  Inch,  0.003  Inch,  and  0.008  Inch 
were  calibrated.  A calibration  curve  for  the  smallest  (0.0011  Inch)  nozzle  Is  shown 
In  Figure  2-6.  This  figure  shows  a region  of  Instability.  Such  regions  are  related 
to  nozzle  diameter  and  flow  rate;  this  relationship  can  be  seen  In  Figure  2-7  for 
the  largest  (0.008  inch  diameter)  nozzle. 

2.3  OPTICS  AND  PHOTOGRAPHY 

In  order  to  satisfy  test  objectives  a photography  system  capable  of  obtaining 
multi  frame  data  over  time  Intervals  on  the  order  of  1 microsecond  was  required.  The 
camera  selected  for  this  study  was  an  Image  converter  having  the  trade  name  Hadland 
Imacon,  manufactured  In  England.  Due  to  late  delivery  of  the  Imacon,  a single  fram- 
ing Beckman  Whitley  camera  was  employed  during  the  early  stages  of  the  prog-am. 

The  Initial  optics  system  is  schematically  shown  In  Figure  2-8.  A spectro- 
physlcs  laser  was  employed  as  a light  source.  The  data  obtained  with  this  system 
was  low  quality.  This  resulted  from  the  fact  that  Interference  due  to  the  coherent 
nature  of  the  source  produced  fringes  on  the  Image  plane;  moreover,  the  optical  sys- 
tem did  not  focus  on  the  droplets,  consequently  any  extraneous  object,  such  as  dirt 
on  the  windows  or  lenses  between  the  light  source  and  Image  plane,  appeared  In  the 
photographs . 

The  system  was  redesigned  for  the  Imacon  as  shown  In  Figure  2-9.  The  Imacon 
Is  electronically  shuttered  and  a multiframing  capability  Is  attained  by  exposing 
sequentially  different  portions  of  a sensitized  tube.  Various  tubes  sensitive  to 
light  ranging  from  Infrared  to  ultraviolet  can  be  employed  with  the  Imacon.  For  the 
present  program  an  S-11  tube  was  used.  Various  plug  In  units  can  be  used  having 
different  framing  rates  and  exposure  times.  The  following  table  lists  the  various 
options  for  recording  eight  16mn  x 18mm  Images  or  sixteen  16mm  x 9 mm  Images. 
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Figure  2-5.  Nozzle  assembly. 
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Figure  2-6.  Droplet  size  - Frequency  relationship. 
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TABLE  2-1.  TRIGGER  MODE  PLUG  IN  UNITS  FOR  HADLAND  IMACON 


Framing  Rate 
(Sec-1) 

Exposure  Time 
(Sec)  x 109 

2 x 106 

100 

5 x 106 

40 

1 x 107 

20 

2 x 107 

10 

The  first  two  units  were  employed  on  the  present  program  giving  total  recording  times 
of  4 usee  (x  10s  frames  per  second)  and  1.6  usee  (5  x 106  frames  per  second). 

A photograph  of  the  test  section,  droplet  generator,  and  Imacon  is  shown  in 
Figure  2-10.  A water  removal  system  was  also  installed  to  eliminate  the  formation 
of  puddles  and  splashing.  The  removal  system  consisted  of  a filter  mounted  directly 
under  the  droplet  stream  flush  with  the  wall  of  the  shock  tube  connected  by  a drain 
line  to  a collection  tank.  The  collection  tank  was  connected  to  a vacuum  pump  and 
instrumented  with  a pressure  gauge.  In  addition  to  eliminating  splashing,  the  water 
removal  system  eliminated  spurious  results,  such  as  were  reported  in  Reference  10, 
arising  from  alteration  of  shock  properties  due  to  the  presence  of  water  vapor. 
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Figure  2-10.  Test  section,  droplet  generator,  camera  assembly. 


SECTION  3 


RESULTS  ANO  ANALYSIS 

Test  results  were  obtained  for  droplet  diameters  ranging  in  size  from  53um 
to  180ym,  velocities  ranging  from  8.2  to  15.5  Kfps,  and  pressure  ranging  from  0.13 
to  0.50  atm.  A summary  of  test  conditions  and  data  for  successful  runs  is  given 
in  Appendix  A.  Droplet  distortion  and  droplet  displacement  data  were  obtained  di- 
rectly from  the  spark  shadowgraphs  obtained  with  the  Hadland  Imacon. 

The  camera  was  triggered  by  the  ion  probe  located  at  station  N with  a time 
delay  based  on  shock  velocity.  This  delay  is  set  prior  to  the  shot  and  is  a func- 
tion of  the  capacitor  voltage  and  pre-shot  shock  tube  pressure.  The  shock  velocity 
is  not  exactly  repeatable  leading  to  errors  in  computed  pre-shot  camera  trigger  de- 
lay times.  This  caused  in  some  cases  early  camera  triggering  resulting  in  frames 
showing  an  unshocked  droplet  and  in  some  cases  late  camera  triggering  resulting  in 
missing  early  time  events.  In  some  cases  for  which  the  delay  was  set  at  too  large 
of  a value  the  droplet  was  swept  out  of  the  field  of  view  by  the  shock  wave  prior 
to  the  completion  of  camera  recording.  It  is  pointed  out  that  difficulties  arising 
from  uncertainty  in  trigger  delay  times  could  be  circumvented  by  increasing  total 
recording  time;  however,  this  results  in  a reduction  of  the  time  resolution  of 
events. 

3.1  EXPERIMENTAL  RESULTS 

Results  based  on  86  frames  and  14  runs  are  tabulated  in  Appendix  A.  Figure 
3-1  is  an  example  of  shadowgraph  data.  From  approximately  80  data  points  of  the 
type  shown  in  Figure  3-1  a sequence  of  droplet  deformation  can  be  constructed;  this 
sequence  is  shown  in  Figure  3-2.  Assigned  to  each  deformation  stage  shown  in 
Figure  3-2  is  a shape  index  varying  from  0 to  7;  0 corresponds  to  the  unshocked 
droplet  and  shape  index;  7 corresponds  to  droplet  fragmentation  or  break-up.  Ini- 
tially a small  amount  of  mass  is  stripped  from  the  droplet  and  it  remains  undeformed 
(shape  index  1).  Subsequently  the  droplet  begins  to  dilate  in  a direction  normal 
to  the  direction  of  motion  of  the  shock  wave  (shape  index  2).  Then  a tail  begins 
to  form  (shape  index  3),  the  tall  widens  and  elongates  (shape  index  4),  and  becomes 
irregular  (shape  index  5);  this  is  followed  by  the  formation  of  waves  on  the  front 
surface  of  the  droplet  (shape  index  6)  and  finally  fragmentation  or  break-up  of  the 
droplet  takes  place  (shape  index  7). 
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Figure  3-1.  Selected  frames  from  three  tests. 


Figure  3-2.  Droplet  deformation. 
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In  order  to  correlate  the  data,  the  following  dimensionless  parameters  are 
Introduced. 

X - x/d0 

t = u2t/d0 

where  x = droplet  displacement 

dQ  * undisturbed  droplet  diameter 

U2  * gas  velocity  behind  shock 

t * time  measured  from  shock  arrival  at  droplet 

P4  * droplet  density 

P2  = gas  density  behind  shock 

Figure  3-3  shows  the  variation  of  the  shape  index  defined  by  Figure  3-2  with 
dimensionless  time  T.  Figure  3-4  shows  dimensionless  displacement  data  plotted  as 
a function  of  dimensionless  time.  The  scatter  in  shape  index  data  is  partially 
attributable  to  the  fact  that  it  is  semi -quantitative  in  nature.  Both  the  shape 
index  data  and  displacement  data  indicate  the  shock  layer  effects  to  be  smaller  at  a 
free  stream  pressure  of  380  torr  than  100  torr;  that  is,  at  380  torr  the  shape  in- 
dex and  displacement  variation  with  dimensionless  time  is  less  than  that  at  100  torr. 
It  is  pointed  out  that  the  present  data  is  in  reasonably  good  agreement  with  the 
correlation  of  Engel’s  data  but  dfffers  significantly  from  AVCO's  correlation. 

The  present  data  indicates,  that  early  time  droplet  displacement  takes  place  at  a 
significantly  lower  rate  than  late  time  displacement.  That  is  during  early  time, 
the  measured  rate  of  displacement  Is  imperceptible,  while  at  later  times  it  varies 
almost  quadratically  with  time.  Moreover,  the  duration  of  this  so-called  early 
time  behavior  is  longer  at  higher  pressures  (380  torr)  than  at  lower  pressures  (100 
torr).  However,  Inconsistencies  existed  between  photographic  and  ion  probe  data  at 
high  pressures;  that  is,  the  ion  probe  data  Indicated  shock  arrival  times  at  the 
droplet  approximately  4 usee  later  than  was  indicated  by  photographic  data.  The 
times  employed  in  Figure  3-3  and  3-4  for  380  torr  data  were  obtained  from  photo- 
graphic rather  than  ion  probe  information.  Due  to  the  inconsistency  mentioned  above 
the  subsequent  analysis  weights  the  100  and  190  torr  data  more  heavily  than  the  380 
torr  data. 

3.2  ANALYSIS 

The  results  presented  in  the  previous  section  of  this  report  have  been  em- 
ployed to  predict  shock  layer  attenuation  of  recession.  In  order  to  do  this  the 
effects  of  shock  layers  on  impacting  kinetic  energy  for  a given  droplet  distribution 
must  be  determined.  The  analysis  is  carried  out  for  conditions  corresponding  to 
those  in  the  stagnation  region  of  a reentry  vehicle. 
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displacement  data. 


In  body  fixed  coordinates  the  Impact  velocity,  u^,  will  be  the  relative  velo- 
city between  the  body  u®  and  the  particle  up  at  Impact  time  t^.  That  Is 

u<  = u®  - un(tj  (3-1) 


The  particle  velocity  as  a function  of  time  can  be  obtained  from  the  derivative  of 
the  displacement  curve  x(t).  Thus 


Up ( t ) = x(t) 


The  dimensionless  impact  velocity  can  be  obtained  from  the  derivative  ot  the 
displacement  curve  shown  in  Figure  3-4  utilizing  the  following  relation 

Vi  - 1 - (1  - X^)  (3-3) 

where  V1  = u^ 

e 1 oJo2 

* density  of  droplet 
U2  - (1  - t)u„ 

and  the  subscripts  1,  °°,  and  2 denote  conditions  at  impact,  in  the  free  stream,  and 
behind  the  shock  wave,  respectively.  For  reasons  mentioned  previously  the  lower 
pressure  data  is  felt  to  be  more  accurate  and  this  data  has  been  weighted  more 
heavily  than  the  380  torr  data  in  the  correlations.  Consequently,  the  following 
correlation  of  displacement  data  has  been  adopted 


0 < T < 0.75 


X « 0.5207  - 1.388  T ♦ 0.9256  V 0.75  < T 


The  above  two  equations  represent  a reasonably  good  fit  to  the  data  while  maintain- 
ing continuity  of  displacement,  X(T),  and  the  dimensionless  time  derivation  of  dis- 
placement, X(T).  As  can  be  seen  from  Equation  (3-3),  the  dimensionless  impact  time, 
Tj,  is  required  In  order  to  obtain  V^.  The  relative  distance  traveled,  0,  between 
the  shock  front  and  the  body  at  time  t is  given  by 

D - u»t  - x(t) 


(3-5) 


and  the  impact  time,  t^,  Is  obtained  by  replacing  D with  the  shock  stand-off 
tance,  A$  giving 

A$  = Ucot.  - 

Consequently  the  dimensionless  impact  time  is  obtainable  from 

A --  T./(l  - c)^7p^  - X(T.) 


where  A = A$/do 

A$  = shock  stand-off 

dQ  = undisturbed  droplet  diameter 

and  the  shock  stand-off  distance  in  terms  of  nose  radius  R^  is  given  by 

A$  = 0.78  e 

Thus  for 


A < .75//p2/ P^l  - e) 

the  dimensionless  Impact  time  is  given  by 

T1  x ‘ 


and  for 


A > - 75//p2/p&( 1 - e) 


„ _ Z - Vz2  - 3 . 702( A + 0.5207) 
Ti ore 


dis- 

(3-6) 
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(3-8) 


(3-9) 


where 


The  present  data  was  obtained  for  conditions  representative  of  a reentry  ve- 
hicle stagnation  region.  For  this  region,  in  contrast  to  that  on  a vehicle  frustum, 
mass  loss  is  dominated  by  catastrophic  breakup  characterized  by  a breakup  time,  Tb, 
rather  than  stripping.  Thus  for  computing  the  impact  of  shock  layers  on  erosion  it 
will  be  assumed  that  the  impacting  mass  is  given  by 

mi  = % 
m^  « 0 

Implicit  in  this  assumption  is  that  the  droplet  fragments  vanish  on  a time  scale  that 
is  small  compared  to  the  time  it  takes  for  the  original  droplet  to  fragment.  This  is 
reasonable  since  the  fragments  are  small  compared  to  the  orginal  droplet  and  since 
they  have  been  heated  for  a significant  period  of  time.  However,  for  a number  of  rea- 
sons this  assumption  is  not  critical.  The  result  of  this  assumption  is  that  the  kin- 

etic energy  of  Impacting  particles  below  a critical  diameter  is  zero.  As  will  be 
shown,  the  kinetic  energy  of  particles  in  the  region  of  breakup  has  already  decreased 
significantly  due  to  drag  effects  and  for  most  flight  conditions  will  become  negligi- 
bly small  for  T^  20  percent  greater  than  Tb  using  the  existing  displacement  correla- 
tion. Moreover,  the  predicted  effects  Integrated  over  a droplet  size  distribution 

will  be  shown  to  not  be  a very  sensitive  function  of  breakup  time  due  to  the  fact  that 

most  of  the  total  mass  entering  the  shock  layer  is  distributed  among  particles  having 
diameters  greater  than  critical  droplet  diameters. 

The  data  shown  in  Figure  3-3  has  been  employed  to  deduce  breakup  times.  Three 
curves  are  shown  plotted  through  the  data,  one  curve  is  a least  squares  fit  through 
the  high  pressure  data,  one  through  the  100  and  180  torr  data  and  one  through  all 
of  the  data.  Associating  breakup  with  a shape  index  of  7 yields  a breakup  time  range 
of  3.5  to  4.0.  For  reasons  mentioned  earlier,  the  high  pressure  data  is  suspect, 
consequently.  It  is  felt  that  the  low  end  of  the  above  range  represents  a more  real- 
istic estimate  of  breakup  time  than  the  high  end.  Thus,  a realistic  estimate  for 
breakup  time  is  3.5  whereas  a conservative  estimate  (minimizing  shock  effects)  is 
about  4.0.  It  will  be  shown  that  the  net  effect  of  the  shock  layer  on  erosion  pre- 
dictions does  not  vary  significantly  over  a range  of  Tfa  from  2.75  to  4.0  for  stagna- 
tion region  conditions,  thus,  accurately  pin-pointing  the  value  of  Tb  Is  unnecessary. 

Predictions  based  on  the  modeling  described  above  are  compared  to  those  of 
AVCO;  hence.  It  Is  appropriate  to  describe  the  methodology  upon  which  the  AVCO  pre- 
dictions are  based.  Impact  velocities  are  obtained  from  the  derivative  of  the  fol- 
lowing correlating  function  (Reference  9 and  12) 

X • 0.8  T2  (3-11) 
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Ti  < Tb  ^ 


T1  >Tb 


(3-10) 


r 


breakup  Is  now  given  by  (Reference  12). 


\ 


TboJ*  - 35/We  ^ (3-12) 

where  We  ■ P2U2do^° 

o * 75  dvnes/cm  (surface  tension  at  room  temperature) 
qm  k = f (p,  u)  = 1 for  present  test  conditions 

and  breakup  Is  determined  by  comparing  from  Equation  (3-7)  to  Tb  from  Equation 
(3-12). 

In  order  to  predict  effects  of  shock  layers  on  recession  rate  the  following 
mass  loss  law  is  employed  with  both  the  AVCO  and  Aerotherm  shock  layer  models 

G * ^uj  (3-13) 

The  recession  rate  is  then  given  by 

s = m(K1u^)/p$  (3-14) 

where  K1  = constant  (material  dependent) 

m « mass  flux  of  Impinging  particles 
p$  * density  of  surface  material 
Uj  ■ Impact  velocity 

Therefore  the  ratio  of  recession  rate  Including  shock  layer  effects  to  that  neglec- 
ting these  effects  (s^/sj  is  given  by  the  following  kinetic  energy  flux  ratio 

il . . a (3-15) 

2.  <«*>.  ’ ’ 

The  above  equations  can  now  be  used  to  obtain  the  (s^/s^)  ratio  as  a function  of 
altitude,  velocity,  nose  radius,  and  particle  diameter.  In  practice  the  droplets 
are  distributed  over  a diameter  range  and  In  order  to  obtain  the  net  effects  of 
shock  layers  In  flight,  the  above  equation  must  be  Integrated  over  some  suitable 
distribution  function.  As  was  mentioned,  It  will  be  assumed  that  all  particles 
with  diameters  below  the  critical  value  d£  do  not  contribute  to  the  erosion  damage 
where  dc  Is  the  diameter  of  particles  having  Impact  time  T^  equal  to  Tb>  The  follow 
Ing  mass  distribution  function  (Reference  24)  Is  used  for  computational  purposes 
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(3-16) 


M(d0) 


where  the  mass  mean  diameter  is  40  . 

m 

The  quantity  (uddQ)  is  the  mass  fraction  of  drops  in  the  diameter  range  be- 
tween dQ  and  dQ  + ddQ.  The  actual  recession  rate  ratio  for  a given  distribution  is 
now  given  by 


• K actual  recession  rate  for  given  droplet  distribution 
recession  rate  without  shock  wave 


l - f «,Vf  u(d0)<W0  - f »»  u,<d0)dd0  (3-17) 

o dc 

The  preceding  equations  have  been  coded  and  predictions  can  be  made  for  arbitrary 
combinations  of  altitude,  velocity,  R^,  Dm,  and  assumed  dimensionless  breakup  time, 

T^.  A few  examples  will  be  presented  for  selected  flight  conditions  and  comparisons 
made  with  predictions  using  AVCO  breakup  times  and  accelerations. 

According  to  the  present  model  the  shock  layer  attenuates  recession  by  (a)  re- 
ducing impact  velocity,  (b)  shattering  the  droplet.  Damage  from  droplets  with  impact 
time,  Tj,  greater  than  breakup  time,  T^,  will  be  reduced  only  due  to  deceleration. 

The  dimensionless  velocity,  V^,  Is  a measure  of  the  deceleration.  It  is  of  interest 
to  compare  dimensionless  velocities  predicted  using  the  AVCO  and  Aerotherm  displace- 
ment correlations.  From  equation  (3-3)  it  can  be  seen  that 


V1  = 1 - F(Alt. , Vel . ) X(Ti) 

(3-18) 

where 

• 

XAVC0  " 1,6  T 

Wo  • °* 

T <0.75 

XAero  " 1>851  T -K388‘ 

T >0.75 
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Table  3-1  gives  values  for  the  quantity  F,  appearing  In  the  above  equation,  at  dif- 
ferent values  of  altitude  and  velocity. 

TABLE  3-1.  VALUES  OF  COEFFICIENT  F APPEARING  IN  EQUATION  (3-18). 


Alt. 
(Kft. ) 

Vel. 

(Kfps) 

F 

10 

10 

0.073 

10 

20 

0.086 

20 

10 

0.062 

20 

20 

0.076 

In  order  to  compare  dimensionless  impact  velocities,  a nominal  value  of  0.074  is 
shown  for  F.  Table  3-2  shows  a comparison  of  impact  velocities. 


TABLE  3-2.  COMPUTED  IMPACT  VELOCITIES  FROM  EQUATION  (3-18) 

FOR  F = 0.074 


T 

Vi 

(AVCO) 

Vi 

(Aero) 

0 

1 

1 

0.5 

0.94 

1 

1.0 

0.88 

0.97 

2.0 

0.76 

0.83 

4.0 

0.53 

0.56 

The  differences  between  Impact  velocity  are  not  large;  however,  their  differ- 
ences will  lead  to  larger  differences  In  recession  rate  due  to  the  fact  that  reces- 
sion rate  varies  as  the  kinetic  energy  of  the  Impacting  droplet. 

The  effect  of  the  difference  In  Impact  velocity  and  break-up  time  on  predicted 
recession  rate  as  a function  of  droplet  diameter  Is  shown  In  Figure  3-5.  The  gradual 
recession  rate  drop  off  with  decreasing  droplet  size  Is  attributed  to  droplet  velocity 
reduction,  while  the  sharp  drop  off  portion  of  the  curve  Is  attributed  to  droplet 
breakup.  The  differences  between  the  present  predictions  and  those  of  AVCO  due  to 
Impact  velocity  are  related  only  to  differences  In  the  displacement  correlation  dis- 
cussed previously.  The  differences  due  to  the  breakup  effect  Is  related  to  differ- 
ences In  two  factors:  (a)  displacement  correlation,  (b)  breakup  time.  This  arises 
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from  the  fact  that  the  criterion  for  breakup  (T\ , > Tb)  involves  both  the  impact  time 
and  the  breakup  time  and  the  impact  time,  , is  computed  from  the  droplet  displace- 
ment correlation.  Impact  times  computed  from  the  AVCO  correlation  are  greater  than 
those  predicted  from  the  present  correlation  due  to  greater  deceleration;  moreover 
AVCO  breakup  times  employed  by  AVCO  are  less  than  present  values.  Hence,  the  two 
factors  mentioned  above,  contributing  to  differences  due  to  breakup  effects,  both 
reinforce  each  other. 

Differences  between  the  AVCO  and  present  values  for  breakup  time  are  attribut- 
able to  differences  in  interpretation  of  shadowgraph  data.  In  the  present  investiga- 
tion breakup  in  assumed  at  times  corresponding  to  the  onset  of  large  waves  (observed 
from  shadowgraphs)  on  the  windward  side  of  the  droplet.  X-ray  absorption  data  obtained 
by  AVCO  apparently  indicates  mass  disappaerance  at  times  prior  to  those  corresponding 
to  the  onset  of  front  surface  waves.  It  is  felt,  however,  that  some  uncertainty  exits 
relative  to  interpretation  of  the  x-ray  data.  It  is  single  frame  data,  thus  requiring 
a number  of  tests  having  identical  test  conditions,  identical  droplet  size,  and  a 
variation  of  times  for  triggering  the  camera  in  order  to  capture  a sequence  event 
leading  to  breakup.  Moreover,  the  shadowgraph  and  x-ray  data  are  not  taken  simulta- 
neously; hence,  in  order  to  relate  the  two,  a process  similar  to  that  mentioned  above 
must  also  be  followed  for  obtaining  shadowgraph  data.  It  is  felt  that  the  possibility 
of  error  arising  from  the  difficulty  in  implementing  both  of  these  processes  in  order 
to  get  one  comparison  between  an  x-ray  picture  and  a shadowgraph  exactly  at  the  break- 
up time  is  significant.  Another  concern  relative  to  the  application  of  x-ray  data  is 
that  it  was  obtained  only  for  relatively  large  droplets  (>1000ym).  In  order  to  ap- 
ply this  data  to  predictions  of  small  droplet  (<100pm)  behavior,  it  is  necessary  to 
have  a high  degree  of  assurance  relative  to  the  validity  of  the  scaling  laws  used. 

Lack  of  understanding  of  the  basic  physics  related  to  breakup  and  a significant  de- 
viation of  data  from  correlation  involving  the  dimensionless  parameters  X,  T,  and 
We  are  indicative  of  the  fact  that  further  work  is  required  relative  to  developing 
scaling  laws.  In  view  of  the  foregoing  discussion  and  for  the  sake  of  conservativism, 
it  is  recommened  that  breakup  times  based  on  the  present  small  droplet  shadowgraph 
data  be  employed  to  predict  small  droplet  breakup  for  flight  conditions. 

Figure  3-6  shows  the  effects  of  nose  radius  on  the  normalized  recession  rate. 

It  can  be  seen  that  for  a 0.5  inch  nose  radius  vehicle  the  shock  layer  offers  little 
protection  while  for  a 2.0  Inch  nose  radius  vehicle  the  protection  is  significant. 

The  results  In  Figure  3-6  have  been  cross  plotted  In  Figure  3-7  to  show  droplet  mass 
mean  diameter  effect  on  net  recession;  this  effect  Is  also  significant.  It  can  be 
seen  that  the  differences  between  the  present  predictions  and  the  AVCO  curves  In- 
crease with  Ru  and  0 . It  It  pointed  out  that  the  results  can  be  presented  In  terms 
of  the  single  parameter  (RN/Dm);  the  previous  plots  were  presented  to  illustrate 
the  effects  separately  of  varying  RN  and  Dm. 
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Figure  3-6.  Effect  of  nose  radius  on 
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Figure  3-8  illustrates  shock  layer  effects  for  relatively  low  altitude  clouds 
with  large  mean  particle  diameter.  Although  the  trends  are  the  same  the  effects  are 

much  less  than  for  the  higher  altitude  case.  Also  for  this  case  the  differences  bet- 

ween the  present  and  the  AVCO  predictions  are  much  smaller  than  for  the  higher  alti- 
tude case. 

A limited  sensitivity  study  was  conducted  in  order  to  assess  the  effect  of 
break-up  time  on  predicted  recession  rate.  The  results  are  shown  in  Figure  3-9. 

The  present  data  indicates  dimensionless  break-up  times  to  be  within  the  range  of 

3-4.  The  sensitivity  of  S to  T^  is  strongly  dependent  on  R^.  For  nose  radii  less 
than  1 inch  the  variation  of  S over  a range  of  Tb  from  3 to  4 is  less  than  6 percent 
and  for  a range  of  T^  from  2 to  3 less  than  19  percent.  For  a nose  radius  of  2 inch 
the  variation  of  S with  T^  is  15  percent  over  the  range  of  T^  from  3 to  4.  These 
results  indicate  that  for  the  present  set  of  flight  conditions,  the  predicted  reces- 
sion rate  is  only  moderately  sensitive  to  break-up  time.  Also  shown  in  this  figure 
are  AVCO  values  of  S for  break-up  times  given  by  Equation  (3-12).  Figure  3-9  illu- 
strates that  AVCO  break-up  times  are  significantly  lower  than  the  present  recommended 
values  of  (3-4).  However,  it  can  also  be  seen  that  for  nose  radii  of  less_than 
1.0  inch  the  difference^  in  break-up  time  do  not  account  for  differences  in  S;  in- 
deed for  = 0.5  Inch  S does  not  vary  with  assumed  T^  at  all.  For  large  RN  the 
predicted  values  do  become  sensitive  to  T^  for  T^  < 3.0.  From  these  results  it  is 
apparent  that  for  nose  radii  less  than  1.0  inch  differences  between  the  present  re- 
sults and  those  of  AVCO  are  attributable  mainly  to  differences  in  acceleration  laws; 
for  larger  nose  radii  the  differences  are  attributable  to  both  differences  in  ac- 
celeration laws  and  differences  in  break-up  time.  It  is  of  interest  to  compare  the 
relative  effects  of  droplet  diameter  on  particle  impact  velocity,  impact  time,  and 
deformation.  These  quantities  are  shown  in  Figure  3-10  and  it  can  be  seen  that  de- 
formation is  significant  over  a much  larger  range  of  dQ  than  break-up  or  accelera- 
tion. For  example  particles  having  diameters  greater  than  60p  will  not  break-up 
and  acceleration  Is  negligible  for  particles  greater  than  120p,  however,  dilation  Is 
significant  (l.e.,  40  percent)  for  particles  as  large  as  300p.  This  observation  is 
also  made  In  Reference  25,  based  In  part  on  Holloman  sled  data. 
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SECTION  4 


CONCLUDING  REMARKS 

The  effects  of  shock  layers  on  droplet  breakup  time  and  acceleration  have  been 
measured.  The  data  has  been  reduced  and  has  been  Incorporated  within  a mathematical 
model  of  the  effects  of  shock  layers  on  stagnation  point  recession  caused  by  hyper- 
velocity flight  through  clouds  of  water  droplets. 

Predictions  of  shock  layer  attenuation  of  recession  rate  were  made  for  selected 
combinations  of  velocity,  altitude,  nose  radius,  and  particle  size.  The  results  were 
compared  to  predictions  made  based  on  a model  Incorporating  AVCO  breakup  time  and 
acceleration  correlations.  The  results  are  summarized  in  the  following  table. 

TABLE  4-1.  SUMMARY  OF  NORMALIZED  RECESSION  RATE  PREDICTIONS 


Alt. 

(Kft) 

Vel . 
(Kfps) 

On.) 

M 

1 

AERO 

AVCO 

12 

13.2 

.5 

20  - 40 

0.77  - 0.92 

0.61  - 0.81 

1.0 

0.43  - 0.77 

0.16  - 0.54 

1.5 

0.20  - 0.59 

0.02  - 0.26 

2.0 

0.08  - 0.43 

0.00  - 0.09 

6 

11.4 

1.0 

100  - 200 

0.94  - 0.97 

0.82  - 0.92 

1.5 

0.87  - 0.95 

0.75  - 0.87 

2.0 

0.79  - 0.94 

0.54  - 0.80 

The  above  table  shows  that  for  the  relatively  high  altitude  small  droplet  size 
case,  shock  layers  begin  to  be  effective  relative  to  reducing  recession  rates  for  1.0 
Inch  nose  radius  vehicles  and  are  extremely  effective  for  2.0  Inch  nose  radius  vehicles. 
For  the  lower  altitude  case  the  effects  shown  are  significantly  smaller;  this  Is  di- 
rectly attributable  to  the  fact  that  larger  mean  droplet  diameters  were  employed  for 
the  low  altitude  cases. 


The  present  model  predicts  smaller  shock  layer  effects  than  does  the  model 
based  on  the  AVCO  data  correlations.  This  Is  attributable  to  (a)  present  droplet 
acceleration  effects  are  smaller  than  those  reported  by  AVCO,  (b)  present  values  of 
breakup  time  are  layer  than  those  reported  by  AVCO.  These  differences  reinforce  each 
other  relative  to  their  Impact  on  recession  rate  predictions.  The  differences  In  ac- 
celeration are  small  (less  than  15  percent);  however,  the  effect  of  these  differences 
on  recession  rate  predictions  Is  significant  (20  to  40  percent)  since  recession  rate 
varies  as  the  kinetic  energy  of  the  Impacting  particle.  The  differences  between  pre- 
sent and  AVCO  values  of  breakup  time  result  from  differences  In  Interpretation  of 
data.  As  previously  described  a more  conservative  approach  was  adopted  In  the  pre- 
sent investigation. 

The  present  results  also  show  that  deformation  is  significant  up  to  much  larger 
droplet  diameters  than  breakup  and  acceleration;  from  this  it  can  be  concluded  that 
droplet  deformation  effects  on  material  response  warrants  investigation. 
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APPENDIX  A 


A summary  of  test  conditions  and  test  results  is  given  In  Table  A-l. 
table  gives  pressure  upstream  of  the  shock,  shock  velocity,  droplet  diameter, 
number,  time  after  passage  of  the  shock  for  each  frame,  droplet  displacement, 
sionless  time,  deformation  shape  index,  and  dimensionless  displacement. 


This 

frame 

dimen- 


TABLE  A.  SUMMARY  OF  TEST  DATA. 


Run 

P 

(atm) 

V 

(kfps) 

r^7_ 

(v) 

Frame 

t 

( usee) 

X 

(mm)xlO3 

T 

S 

X 

26 

0.133 

14.9 

109 

5 

0.02 

0 

0.032 

0 

0 

26 

0.133 

14.9 

109 

6 

0.52 

21.8 

0.851 

2.75 

0.2 

26 

0.133 

14.9 

109 

7 

1.02 

185.3 

1.670 

4.25 

1.7 

26 

0.133 

14.9 

109 

8 

1.52 

436 

2.488 

5.0 

4 

26 

0.133 

14.9 

109 

9 

2.02 

763 

3.306 

6.0 

7 

26 

0.133 

14.9 

139 

5 

0.02 

0 

0.0256 

0 

0 

26 

0.133 

14.9 

139 

6 

0.52 

0 

0.667 

2.5 

0 

26 

0.133 

14.9 

139 

7 

1.02 

55.6 

1.309 

4.00 

0.4 

26 

0.133 

14.9 

139 

8 

1.52 

166.8 

1.95 

4.75 

1.2 

30 

0.122 

15.2 

140 

4 

1.15 

1.44 

3.75 

30 

0.122 

15.2 

140 

5 

1.65 

2.06 

4.5 

30 

0.122 

15.2 

140 

6 

2.15 

2.69 

5.5 

30 

0.122 

15.2 

140 

7 

3.65 

3.50 

7.0 

34 

0.130 

15.6 

180 

2 

0.04 

0 

0.041 

1 

0 

34 

0.130 

15.6 

180 

3 

0.54 

0 

0.557 

2 

0 

34 

0.130 

15.6 

180 

4 

1.04 

45 

1.10 

3 

0.25 

34 

0.130 

15.6 

180 

5 

1.54 

180 

1.59 

4 

1.00 

34 

0.130 

15.6 

180 

6 

2.04 

360 

2.10 

? 

2.00 

36 

0.139 

15.5 

95 

2 

0 

0 

0 

0 

0 

35 

0.139 

15.5 

95 

3 

0.497 

0 

0.996 

2.5 

0 

35 

0.139 

15.5 

95 

4 

0.997 

95 

1.99 

4.0 

1 

35 

0.139 

15.5 

95 

5 

1.497 

380 

3 

6 

4 

48 

0.134 

11.1 

53 

1 

1.20 

1.5 

3.5 

48 

0.134 

11.1 

53 

2 

1.4 

1.75 

4.0 

48 

0.134 

11.1 

53 

3 

1.6 

2.00 

4.5 

48 

0.134 

11.1 

53 

4 

1.8 

2.25 

5.0 

48 

0.134 

11.1 

53 
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